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1. Introduction 


Several recent publications have taken advantage of the improved version of the T-matrix 
method' and presented an e (tensive survey of light-scattering properties of polydisperse, 
randomly oriented spheroids w th moderate aspect ratios (see, e.g., Refs. 2 and 3 and references 
therein). However, the results o F Refs. 4-6 give an indication that spheroids with extreme aspect 
ratios may have scattering prop Arties dramatically different from those of moderately aspherical 
particles. Therefore, in this pc.per we calculate the optical cross sections, the asymmetry 
parameter, and the scattering ma :rix for randomly oriented ice spheroids with aspect ratios up to 
20 and equivalent-sphere size pa ameters up to 12 and compare them with analogous results for 
surface-equivalent spheres and sj heroids with an aspect ratio of 2. Since needle-like and plate- 
like particles can be abundant in various natural and artificial environments, we conclude the 
paper by a discuss on of potential practical implications of our results. 

2. Computations 

We performed calculations for 5 particle shapes: spheres, prolate spheroids with semi-axis ratios 
a/b = 0.5 and 0.05, and oblate spheroids with ajb =2 and 20, where b is the rotation (vertical) 
semi-axis and ci is the horizontal semi-axis of a spheroid. The size of a spheroid was specified in 
terms of the radius of the surface-equivalent sphere r s (or, equivalently, the radius of the sphere 
having the same projected area 7 ). To suppress the interference structure in light-scattering 
patterns (Refs. 2, 3, and 8), the computation results for spheres and spheroids with a/b =0.5 and 
2 were averaged over a narrow gamma distribution of equivalent-sphere radii given by 
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here, his explanation is worth recalling; it requires a brief digrr ision into the theory 
of normal modes. 

A tide of frequency u has a complex tidal admittance Z{u, l ) at location L which 
can be written as an expansion of oceanic normal modes (e.g., Platzman. 1991): 

Z{u,L) = Y l ( l -“/°k)- l S k H k (L). ( 1 ) 

k 

According to this prescription each mode k contributes to tl e tidal admittance an 
amount depending on three factors: a “frequency factor” involving the mode s com- 
plex frequency a k , a “shape factor” depending on the spatial coherence between the 
tidal potential and the mode elevation, and a “location factor” which is simp y the 
modal amplitude H k at L. . 

Normal modes for the world ocean have been computed by several grouts Ac- 
cording to calculations by Platzman et al. (1981) there are three modes within o very 
near the diurnal band, having periods 23.7, 25.7, and 28.7 nours. The latter mode 
is primarily an Antarctic Kelvin wave and Pacific half-wav< , with hardly any i nflu- 
ence in the Atlantic Ocean In contrast, the other two mod^s are primarily A:.k.ntic 
modes. The relevant portion of the 25.7-h mode in shown in Figure 1. It has two 
amphidromes rotating cyclonically in each basin and larges amplitudes in the North 
Atlantic, especially the northeast sector. The 23.7-h mod is not unlike the 25. T -h 
mode in being largest in the North Atlantic and relativel; small elsewhere; it does, 
however, have an additional amphidrome in the equatorial Atlantic. See Platzman et 
al. (1981) and Platzman (1985) for a comprehensive set of mode diagrams. 

As Cartwright (1975) pointed out, the 25.7-h mode is spatially well correlated wi h 
the third-degree component of the diurnal tidal potential, but not with the secon 1- 
degree component. The relevant spherical harmonics of the potential are (Cartwrig it 
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and Tayler, 1971): 


Y 2 l (0, <t>) = 7(5/2470 3 sin 9 cos 9 e i(t> 

Y$ ($,<(>) = 7(7/19270 3 sin0 (5cos 2 9—1) e id> 

for spherical coordinates (d, <f>). Both functions are of order 1, as is appropriate for 
diurnal tides; they are shown in Figure 2. The Y $ harmonic, like the 25.7-h mode, 
is symmetric about the equator, and both the harmonic and the mode have local 
maxima near latitudes 0° and ±60°. In contrast, the Y^ 1 harmonic is anti-symmetric 
about the equator, and it cannot be expected to easily excite a symmetric mode. 
Platzman (1984) confirmed this argument by calculating the correlations: he found 
the shape factor 5* is 0.30 for Y 3 L but only 0.05 for Y^ 1 . (Note that the faster 23.7-h 
mode is also symmetric, but with an additional amphidrome near the equator it is 
less well matched to Y$ than is the 25.7-h mode.) 

If this story of Mi is correct, then one naturally suspects that other third-degree 
diurnal tides are similarly enhanced in the Northeast Atlantic. It is the purpose of 
this note to show that they indeed are. In addition, by studying the third-degree 
tidal admittances across the diurnal band, one may deduce further characteristics of 
the normal modes and their responses, namely their periods — to be compared with 
Platzman 7 s — and their specific dissipation Q ~ l . 

Like Mi the other third-degree diurnal tides are separated from second-degree 
tides by only 1 cycle in 8.8 years. The other third-degree tides are also significantly 
smaller than Mi, which is itself only a few millimeters. Both facts suggest that very 
long time series of hourly data are required before reliable estimates can be obtained. 
The data employed below, from three tide gauges, span at least two complete cycles 
of the lunar node. i.e. at least 37 years. 
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2 Estimates of degree-3 tides 

Hourly sea level data were analyzed from three stations: Newlyn, England (50°6' N, 
5 ° 33 ' W), Lerwick, Scotland (60°6'N, 1°8'W), and Vigo, Spain (42°14'N, 8°44'W). The 
data timespans are: 1960-1997 (Newlyn), 1959-1997 (Lerwick), 1943-1990 (Vigo). 
The data are of high quality and remarkably complete; the longest gaps are a three- 
week period during 1977 at Lerwick and two two-week periods during 1953 at Vigo. 

The tidal analyses were performed by least-squares fits to sinusoids of known fre- 
quencies. All major constituents, including long-period tides and nonlinear compound 
tides, were accounted for, although in the tables given here (Tables 2-4) only the rel- 
evant third-degree diurnals are listed. For these third-degree diurnals, solutions were 
attempted for all lines whosg potential amplitudes exceed 0.3 mm in the Cartwright- 
Tayler expansion. There are five such constituents spanning the diurnal band between 
2Qi J i - 2Qi is more than 10 times smaller than Mi, but it is included because 

of its importance at the edge of the band. The five constituents are separated in 
frequency by one cycle/month, with a single gap where one normally expects the Ki 
constituent (its largest third-degree line is only 0.22 mm, even smaller than 2Qi). 

Because of the extraordinarily long lengths of the time series, relevant nodal side 
lines in constituents can be included as separate estimates, independent of the main 
lines. For the third-degree diurnals, nodal lines above the 0.3 mm cutoff exist for 
constituents Qi, Mi, and J[. The nodal lines of Qi and Ji appear to yield genuine 
third-degree estimates, but the nodal lines of Mi are problematic: one of them clearly 
cannot be separated from second-degree lines, since (see Table 1) lines 155.565 and 
155.645 are of almost identical frequencies, and the other line may be affected by 
nodal modulations in a weak nonlinear interaction between Lo and Ki (equivalent to 
a line at 155.465). Only nodal lines from Qi and J L are included in Tables 2-4. 

The tables give amplitude and Greenwich phase lag and the magnitude of the 
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admittance Z, defined in this context as the ratio of the (complex) elevation to the 
normalized potential. All three quantities include standard error estimates, based on 
the diagonal elements of the leas t-sqYi ares covariance matrix. The admittance errors 
for 2Qi are understandably largest owing to its small amplitude in the astronomical 
potential. The given error estimates appear realistic according to two tests: (1) 
The original time series were partitioned into year-long segments and the statistical 
variability of the yearly estimates was compared with the tabulated standard errors 
(allowing for a y/n factor for n yearly estimates); this test, of course, can only be 
performed for major second-degree constituents, but the agreement in these terms 
appeared satisfactory. (2) The error estimates appear to account adequately for the 
small differences in |Z| and Greenwich phase lag between the Qi and Ji nodal lines 
and their main lines, in keeping with a smooth oceanic admittance. 

The tidal estimates for Mi at Newlyn and Lerwick agree well with those obtained 
by Cartwright (1975) for these same two stations, and station Vigo similarly shows 
an enhanced Mi admittance. In fact, all three stations display a strong enhancement 
across the entire diurnal band. This lends considerable weight to Cartwright’s original 
explanation of a resonance between the tidal potential and a single operating normal 
mode (of Q not too high). Moreover, all five constituents show phase lags increasing 
towards the north, consistent with the normal mode phases seen in Figure L Both 
the enhanced amplitudes and the consistent phases must be counted as yet further 
evidence for the existence of oceanic normal modes, supplementing various other 
evidence discussed by Platzman (1991). 

At all three stations, the admittance at Ji is comparable to or slightly larger 
than that at M[. More interestingly, all three stations show admittance amplitudes 
increasing as the frequency decreases away from Mi, with 2Q t having the largest 
admittance in all cases. If Platzman's single 25.7-h mode is the primary contributor 
to these tides, which is the case according to his tidal synthesis (Platzman 1984, Figure 
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10), then the admittances strongly suggest that his period for this mode is too low. 
A period below that of Q u and possibly below that of 2Qi, appears more consistent 
with the observed tidal admittances.' An error of 10-20% in Platzman’s period is not 
unreasonable given various approximations made in the original calculations (e.g., 
neglect of self-attraction and crustal loading). 

3 Modal fits 

With the tidal admittance defined across the diurnal band, it is straightforward to 
use the formalism of Equation (1) to fit hypothetical normal modes (e.g., Platzman 
1991). In particular, if one dominant normal mode is presumed to account for the 

behavior across the band, then we may fit the admittance Z by a simple equation of 
form 

Z(u>) = (1 - ui/a)~ x b. (2) 

where the modal strength b and frequency a are treated as free parameters. Here 
Z, 6, and o are complex, while w is real. Of primary interest is the frequency a 
which determines the modal period T = 2n/Rea and specific dissipation Q~ l = 

2 Im cr/Re a. 

Figure 3 shows the result of fitting, the tidal admittances at Newlyn, using a 
weighted, nonlinear, least-squares algorithm. In this case, there are 2 complex pa- 
rameters to be determined from 5 complex admittances. The resulting frequency 
corresponds to a modal period of 28.1 hours and Q of 7. This estimate of Q is about 
twice the value determined for the diurnal-band Antarctic Kelvin wave but less than 
the Q determined for the semidiurnal tides of the . North Atlantic (e.g.. Platzman 
1991). It is clear that the fitted period depends sensitively on the admittance at 2Q i: 
since its amplitude falls well short of the linear trend between Oi and Q[, the fitted 
curve is forced downward, causing a peak (and hence estimated period) in the general 
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vicinity of 2Qi. In this region of the diagram, the fit is satisfactory given the sizes of 
the error bars. On the other hand, the amplitude fit at Ji is clearly poor and suggests 
the presence of one or more additional modes. 

Equation (2) has therefore been augmented with a second oscillator. The resulting 
fit to the diurnal admittances is shown in Figure 4. The estimated periods in the case 
are 27.5 and 23.0 hours, somewhat closer to Platzman’s 25.7 and 23.7 hours. The 
corresponding Q , around 10, is higher than the single-mode fit. One should not, 
however, put too much credence to this augmented fit; we are fitting 4 complex 
parameters to only 5 complex admittances, and the very good agreement shown in 
Figure 4 in fact suggests overfitting. Unfortunately, there are no additional tidal 
frequencies that can be used to constrain a two-mode fit. 

Fits to the data at Lerwick and Vigo are quite similar to those at Newlyn. The 
estimated periods and Q values for all stations are summarized in Table 5. I am 
disinclined to make a choice as to the “best” estimated period for the fundamental 
mode, except to say that the available evidence clearly suggests a period longer than 
Platzman’s original 25.7 hours, a period most likely longer than the period of Qj, and 
perhaps in the vicinity of 27 to 29 hours. 


Acknowledgments. The tide gauge data employed here were obtained courtesy of 
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TABLES 


Table 1: Mi tidal amplitudes at Newlyn, England 


Doodson 

number 

Degree 

Frequency 

(°/h) 

Potential* 

(mm) 

Equilibrium** 

(mm) 

Observed*** 

(mm) 

155.445 

2 

14.48520 

1.37 

0.359 


155.455 

2 

14.48741 

7.41 

1.943 


155.545 

3 

14.48985 

0.59 

0.188 


155.555 

3 

14.49205 

3.99 

1.268 

4.8 

155.565 

3 

14.49426 

0.52 

0.165 


155.645 

2 

14.49449 

0.59 

0.155 


155.655 

2 

14.49669 

20.62 

5.407 

1.2 

155.665 

2 

14.49890 

4.14 

1.086 



* Potential normalized by acceleration of gravity, extracted from tables of Cartwright 
and Tayler (1971). 

** Equilibrium amplitude equals potential x y n where 7 „ accounts for elastic body 
tide. 

*** Sea level amplitudes estimated by Cartwright (1975). 
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Table 2: Third-degree diurnal tides at Newlyn, England 


Tide 


2Qi 

Qi 

Oi 

Mi 

Ji 


Doodson 

number 

125.655 
135.545 

135.555 

145.655 

155.555 

175.555 
175.565 


Frequency 

(7h) _ 
12.84964 
13.39181 
13.39402 
13.94768 
14.49205 
15.59008 
15.59229 


Potential* 

(mm) 

CL35 

0.50 

1.28 

0.65 

3.99 

1.46 

0.59 


Amplitude 

(mm) 

0.81 ±0.22 
1.00 ±0.22 
2.82 ±0.22 
0.93 ±0.22 
4.67 ±0.22 
1.84 ±0.22 
0.77 ±0.22 


Phase lag 
(deg) 
205 ± 16 
247 ±12 
249 ±5 
276 ± 14 
275 ±3 
294 ±7 
325 ± 16 


\ 2 \ 


2.31 ±0.62 
1.99 ±0.43 
2.20 ± 0.18 
1.42 ±0.33 
1.17 ± 0.06 
1.26 ± 0.15 
1.30 ±0.37 


* Potential normalized by acceleration of gravity, extracted from tables of Cartwright 
and Tayler (1971). 


Table 3: ThircJ-degree diurnal tides at Lerwick, Scotland 


Tide 

Doodson 

number 

Frequency 

(°/h) 

Potential 

(mm) 

Amplitude 

(mm) 

Phase lag 
(deg) 

\z\ 

2Qi 

125.655 

12.84964 

0.35 

1.35 ±0.18 

266 ±8 

3.84 ±0.51 


135.545 

13.39181 

0.50 

1.66 ±0.18 

307 ±6 

3.31 ±0.35 

Qi 

135.555 

13.39402 

1.28 

4.08 ±0.18 

307 ±3 

3.18 ±0.14 

O! 

145.655 

13.94768 

0.65 

1.66 ±0.18 

337 ±6 

2.54 ±0.27 

Mi 

155.555 

14.49205 

3.99 

7.25 ±0.18 

332 ±2 

1.82 ±0.05 

Ji 

175.555 

15.59008 

1.46 

2.54 ±0.18 

354 ±4 

1.74 ±0.12 


175.565 

15.59229 

0.59 

1.22 ±0.18 

15 ±9 

2.06 ±0.30 


Table 4: Third-degree diurnal tides at Vigo, Spain 


Tide 

Doodson 

number 

Frequency 

(°/h) ' 

Potential 

(mm) 

Amplitude 

(mm) 

Phase lag 
(deg) 

\z\ 

2Qi 

125.655 

12.84964 

0.35 

1.08 ± 0. 15 

179 ±8 

3.09 ±0.42 


135.545 

13.39181 

0.50 

1.02 ±0.15 

229 ±9 

2.04 ±0.30 

Qi 

135.555 

13.39402 

1.28 

3.11 ±0.15 

213 ± 3 

2.43 ±0.12 

o L 

145.655 

13.94768 

0.65 

1.29 ±0.15 

236 ±7 

1.98 ±0.23 

Mi 

155.555 

14.49205 

3.99 

5.48 ±0.15 

232 ±2 

1.37 ±0.04 

Ji 

175.555 

15.59008 

1.46 

1.83 ±0.15 

250 ±5 

1.25 ±0.11 


175.565 

15.59229 

0.59 

0.68 ± 0.15 

264 ± 13 

1.15 ±0.26 
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Table 5: Period and Q from modal fits 


Station 

One-mode fits 


Two- mode fits 


Period 

Q 

Period 

Q 

Period 

Q 

Newlyn 

28.1 

6.9 

27.5 

11.4 

23.0' 

9.2 

Lerwick 

28.2 

6.3 

27.6 

9.2 

23.4 

9.9 

Vigo 

28.7 

6.7 

27.8 

6.9 

24.1 

11.2 


Periods are in hours. 
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FIGURE CAPTIONS 


Figure 1: The 25.7-hour mode in the Atlantic Ocean as calculated by Platzman et 
al. (1981). Solid lines are 60° phase contours with arrowhead indicating direction 
of rotation. Dashed lines are amplitude contours, with shading denoting areas with 
amplitudes exceeding the global rms amplitude. Reproduced from Platzman (1991) 
with permission. 


Figure 2: Contours of the second-degree (top) and third-degree (bottom) spherical 
harmonic components of the diurnal tidal potential. Solid and dotted lines denote 
regions of opposite sign. 


Figure 3. Amplitude (top) and phase (bottom) of the third-degree diurnal tidal ad- 
mittance at Newlyn, England, corresponding to constituent estimates tabulated in 
Table 1. The smooth curve represents a one-oscillator fit to the data. The amplitude 
peaks at the modal period 28.1 hours. The poor fit to Jj presumably indicates the 
influence of one or more additional modes. 


Figure 4: As in Figure 3 but for a two-oscillator fit to the data. 




Figure 1: The 5. 7-hour mode in the Atlantic Ocean as calculated by Platzman et 
aL (1981). Soli l lines are 60° phase contours with arrowhead indicating direction 
of rotation- D&t led lines are amplitude contours, with s had ing denoting areas with 
amplitudes excc ding the global rms amplitude. Reproduced from Platzman (1991) 
with permission 




Figure 2: Contours of the second-degree (top) and third-degree (fc }tt n) spherical 
harmonic components of the diurnal tidal potentiaL Solid and dc ted ices d eno te 
regions of opposite sign. 




Period (h) 


30 29 28 27 26 25 24 23 



Figure 3: Amplitude (top) and phase (bottom) of the third-degree diurnal tidal ad- 
mittance at Newlyn, England, corresponding to constituent estimates tabulated in 
Table 1. The smooth curve represents a one-oscillator fit to the data. The amplitude 
peaks at the modal period 28.1 hours. The poor fit to Ji presumably indicates the 
influence of one or more additional modes. 
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Figure 4: Ab In Figure 3 but for a tw^oscillator fit to the data. 



